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INTRODUCTION 


Many  current  interceptor  systems  use  a  guidance  system  based  on  propor¬ 
tional  navigation.  In  a  proportional  navigation  scheme,  the  direction  of  the 
interceptor  is  controlled  in  such  a  way  as  to  make  the  rate  of  rotation  of  the 
interceptor  velocity  vector  proportional  to  the  rate  of  rotation  of  the  line- 
of-sight  (LOS)  fran  the  interceptor  to  the  target  Cl] ,  [2] .  The  derivation 
of  the  proportional  navigation  equations  is  carried  out  below. 

Fran  time  to  time  in  the  analysis  and  design  of  automatic  control  sys¬ 
tems,  it  is  desirable  or  even  necessary  to  have  a  closed  form  analytic  solu¬ 
tion  for  nominal  system  response.  Determination  of  optimal  control  laws  and 
stability  studies  often  require  that  such  nominal  solutions  be  available.  In 
the  present  paper,  we  obtain  a  closed  form  approximate  solution  to  the  classi¬ 
cal  proportional  navigation  problem.  This  solution  is  obtained  through  the 
application  of  quasilinearization  to  the  governing  nonlinear  differential 
equations.  As  it  will  be  demonstrated,  this  analytic  solution  exhibits  very 
good  accuracy  when  compared  to  nunerical  solutions  of  the  nonlinear  equations 
of  motion. 

DERIVATION  OF  EQUATIONS 

In  the  following  paragraphs,  we  derive  the  differential  equations  govern¬ 
ing  the  trajectory  of  an  interceptor  following  a  proportional  navigation  gui¬ 
dance  law.  The  geometry  of  the  target  and  interceptor  are  shown  in  Figure  1. 
The  interception  is  assumed  to  take  place  in  a  plane.  The  position  of  the 
interceptor  relative  to  a  fixed  reference  frame  is  given  by  the  position  vec- 


tor  R.  The  target  is  located  at  F^,.  The  position  of  the  target  relative  to 
the  interceptor  is  r.  The  velocities  of  the  interceptor  and  target  are  V  and 
VT  respectively.  As  shown  in  Figure  1.,  the  heading  angle  of  the  interceptor 
is  i  and  that  of  the  target  is  4^,.  The  angle  between  the  velocity  V  and  the 
LOS,  r  is  4  .  Since 


r 


R  , 


(1) 


(2) 


Writing  the  relative  velocity  vector  v  in  polar  coordinates,  we  find: 


v 


+  r($  +  f)e 


e 


(3) 


where,  e  and  ea  are  radial  and  transverse  unit  vectors.  By  inspection  of 
r  o 

Figure  1, 


v  =  [VTCos(4  +  *  -  *T)  -  VCos*]er  -  ty^inf*  +  ♦  -  *y)  -  VSin*]ee  . 


(4) 


Comparison  of  Equations  (3)  and  (4)  shews  that 


r  =  VTCos(*  + 


*  - 


♦T)-VCos* 


(5) 


and 


♦  +  +  =  -V^inU  +  *  -  ♦.p)  +  VSin* 


(6) 


Proportional  navigation  is  defined  by  the  choice  of  the  control  variable 

t  as: 

♦  -  Xf  (7) 


In  this  equation  X  is  called  the  navigation  constant.  Clearly  if  X  is  a  con¬ 
stant,  Equation  (7)  admits  a  first  integral  of  the  form 

*  =  *q  +  *(♦  "  (8) 


where  the  subscript  o  indicates  that  the  variables  are  evaluated  at  sane  ini¬ 
tial  instant.  Fran  Equation  (8)  we  can  write 
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♦  +  ♦  *  a*  +  0 


where  a  =  l  +  Aand$  =  *o-  A^.  l£on  substitution  of  Equations  (7),  (8)  and 
(9)  in  Equations  (5)  and  .(6)  we  obtain: 


r  =  VTCos(o*  +  8  -  *T)  -  VCos* 


r(l  +  A)*  =  -V^inU*  +  0  -  *  )  +  VSin* 


For  the  purposes  of  this  paper,  we  transform  the  two  above  equations  to 
make  r  the  independent  variable.  This  is  easily  accomplished  by  dividing 
Equation  (11)  by  Equation  (10)  and  by  taking  the  reciprocal  of  Equation  (10). 


dt  _ _ 1 

dr  Tv^CosUt  +  0  -  *  )-VCos*»] 


i 


d*  1  +  3  “  *T)  +  VSin*] 

dr  “  r(l  +  X)[VTCos(a*  +  $  -  -  VCos+1 


rcRivATiai  of  algorithm 


In  the  next  few  paragraphs,  an  algorithm  for  obtaining  an  approximate 
analytic  solution  to  the  equations  derived  in  the  previous  section  will  be 
developed.  'This  algorithm  is  based  on  a  technique  called  quasilinearization. 
Ouasil inear izat ion  is  an  iterative  technique  which  can  be  applied  to  the  solu¬ 
tion  of  nonlinear  differential  equations.  The  technique  was  apparently  intro¬ 
duced  and  developed  by  Bellman  and  Kalaba  [3] .  Using  this  technique,  the 
solution  of  a  non-linear  Two  Point  Boundary  Value  Problem  (TPBVP)  is  reduced 
to  the  recursive  solution  of  a  set  of  linear  TPBVP°s.  The  convergence  of  this 
technique  can  be  shown  to  be  quadratic.  This  is  in  contrast  to  the  linear  con¬ 
vergence  rate  of  the  well  known  Picard  iteration. 

In  that  which  follows,  the  quasilinearization  algorithm  will  be  briefly 
derived.  For  further  explanation  of  the  ideas  utilized,  the  reader  is 
directed  to  the  work  of  Bellman  and  Kalaba  previously  cited.  While  the  field 
of  quasilinearization  is  quite  broad  and  rich  in  both  applications  and  theory, 
the  central  idea  lies  in  the  approximation  of  certain  non-linear  functions  by 
means  of  Taylor  series  expansions  in  function  space.  It  should  be  noted  that 
this  type  of  expansion  is  not  unique  to  the  work  of  Bellman,  et  al.  However, 
it  seems  clear  that  the  technique  utilized  in  this  paper  should  be  attributed 
to  Bellman.  * 


We  shall  restrict  the  present  solution  to  cases  where  v  and  are  speci¬ 
fied  constants  and  $  =  0.  These  restrictions  are  consistent  with  approxima¬ 
tions  usually  made  in  the  analysis  of  control  systems.  The  second  assumption 
simply  reflects  a  rotation  of  the  reference  coordinate  system.  In  cases  where 
the  heading  angle  of  the  target  is  constant,  there  is  no  loss  of  generality  in 
utilizing  this  assumption.  Introducing  these  assumptions  into  Equations  (13) 
and  (12),  we  write 


r(l  +  X)|*  *  g(*) 


(14) 


and 


dt 

dr 


d(*) 


(15) 


where 


-VrnSin(a+  +  P)  +  VSinf 
=  VTCos(a*  +  3)  -  VCos* 


(16) 


and 

S  VTCos(af  +6)-  VCos* 


(17) 
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Note  that  if  the  solution  for  ^(r)  can  be  found,  the  solution  for  t=t(r)  can 
be  found  by  quadrature  of  Equation  (15). 

We  begin  the  derivation  of  the  solution  by  expanding  the  right  side  of 
liquation  (14)  in  a  Taylor  series  in  function  space  ,  i.e.. 


r(l  +  X 


)... 


(18) 


where  we  have  truncated  after  the  linear  term.  In  this  equation,  the  sub¬ 
scripts  are  used  to  indicate  the  iteration  number.  For  simplicity  let 


h 


=  13 

"  3* 


By  differentiation  of  Equation  (16)  we  easily  find: 


h  = 


-.v?  -  V2  +  (l+a)WTCos(X*  +  e) 
[VTCos(a^  +  e)-VCos*]2 


(19) 


Utilizing  the  definition  of  h,  we  can  write 


8 


r 

Equation  (18)  as 


r(1  +  X)dr2  "  hn-l*n  =  gn-l  “  VlVl 


(20) 


With  this  equation  in  mind,  the  recursive  algorithm  becomes  apparent.  That 

is,  a  solution  *  is  assumed  which  satisfies  the  terminal  condition 
To 

2 

♦(Rg)  =  Tg.  This  solution  is  used  to  ocnpute  the  functions  gg(r)  =  g(t  =  tg) 
and  hQ(r)  =  h(^  =  +Q).  Those  explicit  functions  of  r  are  then  substituted 
into  Equation  (20)  and  the  differential  equation  solved  subject  to  the  initial 
condition  ^(Rg)  =  Tg.  This  process  is  repeated  recursively  until  satisfac¬ 
tory  convergence  is  obtained. 

As  previously  noted,  if  an  analytic  solution  for  ^  is  known,  the  solu¬ 
tion  for  t(r)  can  be  obtained  by  quadrature  of  Equation  (15).  Note,  however, 
that  the  functional  form  of  the  right  hand  side  makes  it  difficult,  if  not 
inpossible,  to  obtain  the  integral.  Again,  recourse  to  the  technique  of 
quasilinearization  provides  a  recursive  algorithm  for  obtaining  the  desired 
solution.  We  begin  by  expanding  Equation  (15)  in  a  Taylor  series  in  function 
space. 


dt 
_ n 

dr 


d(* 


n-1) 


3d 

3* 


*  (+n  " 

Vl  n 


pn-l) 


2.  Note,  Tg  is  the  LOS  angle  at  the  initial  position  r  *  R. 


(21) 


9 


As  before,  the  subscripts  indicate  the  iteration  nunber.  Again,  for  simpli¬ 
city  let 


e(+) 


_  3d 

=  a*>  * 


By  differentiation  we  obtain: 


a\L£in(af  +  fi)-VSin*p 

e(*)  = - - - 5  .  (22) 

(V^Cos  ( a  f  +  e)~VCos*] 


Thus,  we  are  able  to  write 


dt 
_ n 

dr 


dn-l  +  en-l(*n  " 


'n-1 


)  + 


(23) 


The  above  equation  defines  a  simplified  quadrature  algorithm  for  the  determi¬ 
nation  of  the  distance-to-go  and  time-to-go.  This  equation  is  integrated  sub¬ 
ject  to  the  terminal  condition  t(0)  =  0. 


sorancN 


Having  derived  the  algorithm  in  the  previous  section,  we  now  turn  to  the 
application  of  this  algorithm.  In  this  section  we  shall  assume  an  initial 
solution  (subscript  0)  and  obtain  the  next  approximation  (subscript  1).  Hie 
accuracy  of  the  solution  will  be  discussed  in  the  following  section. 

Hie  quasilinearization  algorithm  requires  that  an  initial  approximation 
to  the  solution  be  assumed.  Obviously,  it  is  advantageous  that  the  initial 
approximation  be  as  accurate  as  possible.  Our  choice  for  the  initial  approxi¬ 
mation  is  based  on  the  well  known  constant  bearing  solution  [1] . 

It  is  easily  shown  that  for  a  constant  speed  target  moving  on  a  straight 

line  path  ,  it  is  possible  to  choose  a  straight  line  interceptor  ttajectory 

2 

which  will  intercept  the  target  .  Hiis  type  of  interception  is  known  as  a 
constant  bearing  interception  since  in  such  cases,  the  angle  between  the 
interceptor  velocity  vector  and  the  LOS  remains  constant.  Note  that  if  the 
initial  engagement  conditions  in  a  proportional  navigation  interception  are 
such  that  the  angular  rate  of  the  LOS  is  zero,  a  constant  bearing  trajectory 
will  result. 

For  our  purposes,  we  choose  the  initial  approximation  +Q(r)  to  be  given 
by 


+0(r)  =  *0 


(24) 


2.  Providing  that  the  interceptor  speed  is  greater  than  that  of 
the  target. 


11 


where,  of  course,  is  the  initial  bearing  angle  to  the  target.  When  Equa¬ 
tion  (24)  is  substituted  into  Equations  (16),  (17),  (19),  and  (22),  the  fol¬ 
lowing  constants  are  defined! 


^  “Vin(oVe)+vsirtfo 

^  V  5  Go  =  vTcos(o'r0+e)-vcos¥£)  (25) 


W  s  Ho 


-a^-V2+(l+a)WTCos(  A¥0+B ) 
[VTCos  ( at  0+B  J-VCosfp] 2 


(26) 


d0(V  " 


VTCos(a¥0+B)-VCos¥0 


(27) 


aV^5in(a¥0+B )-VSirfr0 
(VTCos ( a*0+B )-VCos¥0) 2 


(28) 
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When  the  above  four  constants  are  substituted  into  Equations  (20)  and 
(23)  we  obtain  the  following  constant  coefficient,  linear  differential  equa¬ 
tions: 


r(ln)  -  Vx  -  G0  -  H0Y0 


(29) 


and 


3T  -  po  +  Eo(*rV 


(30) 


When  these  two  equations  are  solved  subject  to  the  initial  conditions, 
T(Fq)  =  Tq  and  t(Fp)=0  ,  we  obtain 


H_ 


w^,I(^,1+A  - 11 


(31) 
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and 


V*b  lDo 


E0G0 


1  (^7  -  1) 
*0 


(32) 


,E0G0 

^hT 


H0+o 


(Hq+o) 


Kjf> 

*0 


-1] 


Equation  (32)  can  be  used  to  obtain  the  time-of-f  light  by  substituting  r 
-  0.  Thus 


T  =p  [5^0  -  D  ] 

Tf  R0  H„+a  D0J 


(33) 


RESULTS 

In  order  to  gain  insight  into  the  accuracy  of  the  first  iteration  solu¬ 
tion,  a  series  of  ocrparisons  of  predictions  made  by  the  quasilinearization 
solution  and  numerical  integration  of  the  exact  equations  has  been  carried 
out.  Some  of  these  ocrparisons  are  reproduced  in  this  paper.  It  should  be 
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noted  that  the  number  of  parameters  in  the  problem  makes  it  difficult  to  draw 
general  conclusions.  Instead,  we  present  specific  results  and  attempt  to 
infer  more  general  results. 


Table  1  shows  a  comparison  of  a  trajectory  computation  using  the  exact 
solution  and  the  guasilinearization  solution.  To  obtain  this  solution,  the 
following  parameters  were  used: 

V  =  900  m/sec 
VT=  300  m/sec 
A  =4 


Note  that  the  exact  equation  for  $  has  a  singularity  at  r=0.  Thus,  the  numer¬ 
ical  solution  is  not  reliable  in  the  vicinity  of  r=0.  For  this  reason,  the 
last  entry  column  3  of  Table  1  is  indicated  by  asterisks. 


TABLE  1.  Comparison  of  Exact  and  Quasilinearization  Solutions  of 
an  Interception  Trajectory 

PSI  t 


r 

QUASI 

EXACT 

QUASI 

EXACT 

1.000 

-.0146 

-.0144 

.0000 

.0000 

.900 

-.1523 

-.1523 

1.6401 

1.6401 

.800 

-.3195 

-.3194 

3.2831 

3.2830 

.700 

-.5052 

-.5049 

4.9293 

4.9288 

.600 

-.7145 

-.7137 

6.5791 

6.5777 

.500 

-.9554 

-.9535 

8.2330 

8.2298 

.400 

-1.2407 

-1.2365 

9.8917 

9.8854 

.300 

-1.5933 

-1.5845 

11.5561 

11.5448 

.200 

-2.0630 

-2.0439 

13.2279 

13.2081 

.100 

-.7971 

-2.7503 

14.9102 

14.8755 

0.000 

-8.9066 

16.6155 

16.5449 

Examination  of  Table  1  shows  that,  for  this  case  at  least,  the  quasil¬ 
inearization  solution  provides  an  accurate  solution.  In  general,  if  the  ini¬ 
tial  LOS  angle  was  small,  guasilinearization  provided  an  accurate  solution. 
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As  the  initial  LOS  angle  became  much  greater  them  the  LOS  angle  required  for  a 
constant  bearing  solution,  the  accuracy  of  quasilinearization  solutions 
degraded  somewhat. 

Clearly,  the  error  in  the  solution  is  greatest  at  the  end  of  the  trajec- 
tion  (r=0).  It  would  thus  appear  that  a  reasonable  measure  of  accuracy  would 
be  the  percent  error  of  the  time  of  intercept ion.  A  number  of  solutions  were 
carried  out  to  obtain  the  set  of  error  curves  shown  in  Figure  2.  In  this  fig¬ 
ure,  the  independent  variable  is  the  initial  LOS  angle  (To).  The  initial 
azmiuth  angle  (measured  clockwise  from  a  direction  parallel  to  the  target 
velocity  vector)  was  varied  parametrically  to  generate  the  family  of  plots. 
Examination  of  this  figure  shows  that  the  maximum  error  in  the  use  of  the 
quasilinearization  solution  remains  remarkably  small  over  a  substantial  range 
of  trajectory  conditions.  It  appears  that  for  practical  cases  where  engage¬ 
ments  would  be  made  with  LOS  angles  within  ten  or  fifteen  degrees  of  a  con¬ 
stant  bearing  solution,  the  error  induced  by  using  the  approximate  solution 
will  be  less  than  two  or  three  percent. 

CONCLUSIONS 

Having  derived  an  approximate  solution  for  a  proportional  navigation  tra¬ 
jectory  in  an  earlier  section  and  having  corpared  the  approximate  solution  to 
ntmerical  integrations  of  the  exact  equations  in  the  last  section,  we  are  in  a 
position  to  ocntnent  on  the  efficacy  of  the  solution  technique.  As  was  previ¬ 
ously  noted,  in  cases  where  the  initial  LOS  angle  is  relatively  near  the 
lead  angle  for  a  constant  bearing  solution,  the  quasilinearization  solution  is 


very  accurate.  As  this  angle  becomes  larger,  the  accuracy  degrades. Since  the 
quasilinearization  procedure  produces  a  sequence  of  solutions  which  converges 
quadratically,  an  additional  iteration  will  provide  even  more  accuracy. 

Since  in  certain  aspects  of  control  system  design,  it  is  advantageous  to 
have  an  explicit  approximate  solution,  the  use  of  quasilinearization  can  be  of 
real  benefit.  In  such  control  applications,  it  is  usually  desired  to  use  time 
as  the  independent  variable.  Ihe  present  solution  utilizes  range-to-go  as  the 
independent  variable.  Since,  in  general,  the  first  iteration  results  in  a 
nonlinear  function  of  range,  analytic  inversion  of  the  solution  to  functions 
of  time  is  not  practical.  However,  the  form  of  the  solutions  is  simple  enough 
that  numerical  inversion  may  be  practical  even  in  on-line  applications. 

Overall,  it  is  clear  that  the  use  of  the  first  iteration  of  the  quasil¬ 
inearization  solution  for  proportional  navigation  provides  a  straightforward 
and  accurate  approximate  solution.  In  a  subsequent  paper,  we  shall  explore 
the  relative  accuracy  of  the  second  iteration. 


REFERENCES 


[1]  Besserer,  C.  "MISSILE  ENGINEERING  HANDBOOK,  " 

Van  fJostrand,  Prinoeton,  N.J. ,  1958. 

[2]  ns  Army,  "ORDNANCE  ENGINEERING  DESIGN  HANDBOOK  SURFACE  TO  AIR  MISSILE 
SERIES  PART  TWO,  "  OREP  20-292,  Jun  62. 

[3]  Bellman, R.  and  Kalaba,  R.  "CUASILINEARIZATTOH  AND  NONLINEAR  BOUNDARY 
\WJ)E  PROBLEMS,  "  Anerican  Elsevier  Publishing  Co.  Inc,  Hew  York,  1965. 


19 


DISTRIBUTION 


Dr.  Dale  W.  Alspaugh 
Purdue  University 
North  Central  Campus 
Westville ,  IN  46391 

IIT  Research  Institute 
ATTN :  GACIAC 

10  West  35th  Street 
Chicago,  IL  60616 

DRSMI-LP,  Mr.  Voigt 

DRSMI-RDF 

DRSMI-RPR 

DRSMI-RPT  (Record  Copy) 
DRSMI-RPT  (Reference  Copy) 


